Evolving Stellar Background Radiation 
and Gamma-Ray Optical Depth 
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Abstract. We present a semi-empirical model for the evolving far- infrared to ultra- 
violet diffuse background produced by stars in galaxies. The model is designed to 
reproduce the results of deep galaxy surveys, and therefore may be considered as a 
cosmology-independent lower limit to the extragalactic background light. Using this 
model and recent HEGRA data, we infer the intrinsic spectrum at multi-TeV gamma- 
ray energies for Mkn 501 and find that it is consistent with a power law of spectral 
index 2.49 ± 0.04. In turn, this finding renders it rather unlikely that the present-day 
infrared background has an intensity as high as claimed by Finkbeiner et al. [1] . Future 
10 GeV to TeV observations could be used to either constrain the ultraviolet-to-infrared 
background model at high redshifts or cosmological parameters. 



INTRODUCTION 



High-energy gamma-rays originating from sources at cosmological distances can 
be absorbed by pair production in collisions with low-energy photons. Most of the 
low-energy photons are produced by stars in ordinary galaxies and contribute to the 
extragalactic background light (EBL). The role of other contributers, such as Active 
Galactic Nuclei, is under dispute, and can safely be assumed to be unimportant in 
terms of the global energy budget. To model the stellar component of the EBL, we 
use an approach similar to that used in Refs. [2-4], however with the star formation 
rate (SFR) inferred from deep multi-color galaxy surveys. Other approaches aim at 
finding an ab initio description of the SFR from the theory of structure formation 
[5]. Knowing the EBL as a function of redshift, one can calculate the optical depth 
with respect to pair creation and determine the propagation length of gamma-rays 
from cosmologically distributed sources such as Blazars or Gamma Ray Bursts [2,5]. 




FIGURE 1. Left panel: Evolution of EBL from < z < 5 (Az = 0.2). (£:i3^,=0.14, 
ri=150 K, T2=35 K, a=3.5, (3=-l, Zp=1.3, p^{zp)=0.15). Right panel: Present-day EBL with 
(solid line) and without [long dashed line) absorption. Data are Refs. from [10] (open triangles 
and squares /RST), [ll]{circles) , [12]{filled circle at 3.5 • lO^A/DIRBE), [13]("x"), [U]{triangles) , 
[l](squares), [15\{diamond /IRAS, filled circles /DIRBE), [l6\{open circles), [17]{shaded band, EZ- 
RAS detection 1.25 • 104A-50 • lO^A). 



EBL-MODEL 

We start with the spectrum of an evolving simple stellar population (SSP) [6] 
(solar metalicity, Salpeter IMF 0.1-100 Mq). Next, we include absorption and re- 
emission effects (dust and gas of ISM) using a SMC-like extinction law {Eb-v) and 
two modified blackbody spectra for the re-emission from dust (Ti > T2). The Star 
Formation Rate (SFR) is p*(z) oc (1 -|- 2;)"'^ with slopes a for 2; < Zp and /? for 
z > Zp, and a peak value p*{zp). By choosing proper parameters, the emissivities 
ei,{z) obtained from integrating the SFR- weighted SSPs are brought to agreement 
with the data. A second integration from z to 2;max yields the EBL comoving power 
spectrum (Fig. 1) which is independent of cosmology due to a cancellation of the 
cosmology-dependent terms. The EBL increases from high redshifts towards the 
present, but saturates at a roughly constant level beyond the peak in the SFR. The 
EBL spectrum at high z is dominated by young, UV-emitting massive stars. 70% 
of the energy released by stars is re-emitted by dust in the infrared and far-infrared 
bands. The spectrum changes towards its present shape as a result of a growing 
admixture of aged stars. About 40% of the integrated present-day EBL is contained 
in the IR (between 10^ A and 10^A)(see Fig.l and for details Ref. [7]). 



GAMMA-RAY ATTENUATION 



The pair creation optical depth t^^{E, z) for gamma-rays at cosmological redshift 
z and observed energy E is defined as the ratio of coordinate distance and mean 
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FIGURE 2. left panel: Pair-creation optical depth for gamma-rays for a range of redshifts 
(0.03<z<5). right panel: Gamma ray horizon defined by t{z,E) = 1 for various models. 

free path [2]. The optical depth evaluated for our EBL and a range of redshifts 
(0.03 < z < 4) is shown in Fig.2 {Hq = 50 km s"^ Mpc'^, fl = 1, and = 0). 

We can use the optical depth of gamma-rays to infer the intrinsic spectrum of 
MknSOl from recent data. This is achieved by simply multiplying the data with 
the factor exp [xy^ (0.034, i?)], and the result is shown in Fig. 3. 

The gamma-ray horizon t^^{z,E) = 1 defines an interesting relation between 
distance and cut-off energy. In contrast to the EBL, the gamma-ray horizon de- 
pends on cosmology through the integration over a cosmic line element dl/dz. If 
we switch from a standard Cold Dark Matter model (Fig.2, right panels solid line) 
to one with a lower matter density (dashed line, I^a = 0.7), the horizon moves 
closer and the number of sources which can be detected above a given gamma-ray 
threshold energy decreases. The effect could also be mimicked by a flatter SFR [8] . 

CONCLUSIONS 

We have developed a simple model for the evolution of the cxtragalactic back- 
ground radiation field using a minimum number of parameters by fitting the cmis- 
sivity predicted by the model to recent results of high-redshift galaxy surveys. Fu- 
ture improvements will include metalicity effects and a PAH feature around 10//m. 
The intrinsic spectrum of Mkn501 inferred from our model is in good agreement 
with theoretical predictions (see Ref. [9] for a discussion) rendering claims for an 
extremely high level of the infrared background light intensity implausible. The 
EBL model is a powerful tool to analyze gamma-ray spectra from high-redshift 
sources in the 10 GeV to TeV regime for which the pair attenuation occurs in inter- 
actions with ultraviolet-to-infrared photons allowing to probe the era of maximum 
star formation or cosmology. 
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FIGURE 3. Intrinsic spectrum of Mkn501 inferred from our EBL model [upper stars) and fitted 
by a power law with index 2.49 ± 0.04 {upper line). HEGRA original data from Ref. [18]{lower 
stars) and power law times exponential fit {lower line) from Ref. [19] are shown for comparison. 
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